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The  Solution  of  Some  Turning  Point  Problems 


(Abstract) 

This  report  solves  some  turning  point  problems  for  certain 
systems  of  the  form 

^  u  =  [pA(x,p)  +  B(x,p)]u 

where  p  is  a  large  parameter  and  A(x,p)  and  B(x,p)  are  square 
matrices  analytic  in  x  and  p  at  x  =  0  and  p  =  co  . 

The  basic  approach  is  to  obtain  connection  formula  for 
u  in  terms  of  solutions  to  the  comparison  equation 

^  G  =  pA(x,p)G 

For  this  purpose,  define 

^  E  =  pJE 
dx     ^ 

and  , 

V  =  GE'-^ 

and  let 

u  =  Vw 


Then 


^  w  =  [pJ  +  V"^BV]w 


Assume  that  J  is  diagonal  and  in  some  sector  (of  the  complex 
x-plane )  with  vertex  at  x  =  0  the  elements  of  J  and  xV~  BV  are 
HBlder  continuous.   Connection  formula  for  w  in  the  sector  are 
then  obtained  by  an  elaboration  of  a  scheme  due  to  K.  0. 
Prledrichs . 

Two  natural  comparison  systems  are 


^.  /u 
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T-  G  =  p 
dx     '^ 


and 


^^=P 


a. 
a 


/n) 


V  G 


i  =  n 


i  /n 


\ 


G 


1  =  n 


where  each  a.,  is  a  non-negative  integer.   In  each  case  the 
elements  of  G  can  be  written  explicitly  in  terms  of  generalized 
hypergeometric  functions.   These  comparison  systems  are  applied, 
respectively,  to  a  second  order  equation  and  a  fourth  order 
equation  of  the  type  which  arises  in  the  theory  of  hydrodynamic 
stability.   In  each  case  the  integrability  condition  on  V"  BV  is 
verified.   A  typical  result  is  to  obtain  connection  formula  for 
the  second  order  equation 


dx' 


'J 

u  =  [p  a(x)  +  pb(x)  +  c(x,p)]u 


where  the  order  of  the  zero  of  a(x)  and  b(x)  at  x  =  0  is, 
respectively,  a  and  p,  in  terms  of  connection  formula  for 


dx' 


^2  a 
■p  g  =  P  X  g 


whenever 


a 


P  >  ^  - 


I.  Introduction 

The  solution  of  turning  point  problems  via  the  connection 
approach  (which  will  be  described  in  a  moment)  proceeds  in  two 
parts.   (See,  for  example",  V/asow  [2,Chap.6]). 

In  the  first  part  one  obtains  connection  formula  for  solu- 
tions of  a  special  system  by  a  contour  integral  approach  (Watson 
[l.Chap.?]),  or  by  a  Laplace  transform  approach  (Wasow  [1], 
Llnstone  [1],  Lin  [1])  or  by  a  Taylor  series  approach  (Scheffe 
[1],  Turrittin  [1],  Langer  [2]). 

In  the  second  part  first  the  special  system  is  converted 
via  a  1-1  co-ordinate  transformation  so  as  to  bear  a  "suffi- 
ciently close"  resemblance  to  the  original  system.   The  system 
which  results  from  this  transformation  is  called  a  comparison 
system.   Finally,  with  the   aid  of  a  comparison  technique, 
connection  formula  for  the  original  system  are  obtained  in  terms 
of  connection  formula  for  the  comparison  system. 

Needless  to  say,  the  finer  the  comparison  technique  the 
cruder  can  be  the  comparison  system.   In  this  direction,  an 
important  improvement  in  the  standard  comparison  technique  was 
indicated  by  Priedrichs  [l,§22]  where,  for  simplicity,  he  con- 
sidered a  system  of  two  first  order  equations.   An  elaboration 
and  exploitation  of  this  improvement  forms  the  main  theme  of 
this  report.   However,  before  we  discuss  the  improved  comparison 
technique  let  us  first  briefly  review  some  of  the  vital  facts 
concerning  turning  point  problems. 

This  refers  to  Chapter  6  in  the  second  publication  listed 
under  fia.so'M   in  the  Bibliography  at  the  end  of  this  report. 


The  basic  system  to  be  considered  is 

(1.1)  ^  ^  =  ipA(x,p)  +  B(x,p)]u 

where 

u  =  a  column  vector  of  dimension  n 

Ipl  =  a  large  parameter 

X  =  the  complex   independent  variable  ranging  over  a 

"iemann  surface  (R. 

A(x,p)  ,B(x,p)  =  square  nxn  matrices  analytic  in  x 

and  p  at  X  =  0  and  p  =  oo  . 

The  basic  question  to  be  answered  is,  "How  can  one  deter- 
mine the  asymptotic  behavior  (|p|  — >  oo  )  of  a  typical  solution 
of  (1.1)  for  all  X  in//^?"   Complications  arise  because  of  the 
occurrence  of  so-called  turning  points  (sometimes  called  transi- 
tion points)  and  in  order  to  grapple  with  the  heart  of  the  prob- 
lem region  <^  will  he  taken  to  have  a  turning  point  at  x  =  0  and 
no  turning  point  for  x  ^  0,      What  this  implies  for  (1.1)  can  best 
be  described  in  terms  of  the  Jordan  form  of  the  leading  matrix 
A(x,oo).   This  Jordan  form  will  be  assumed  to  be  diagonal  for  x 
in  ^  but  X  /  0. 

Let  \.(x),  1  =  l,2,...,n  denote  the  eigenvalues  of  A(x,oo) 
for  X  ^  0    and  define 

(1.2)  J(x)  =  \S^j.A^(x); 

where 

(O      ,  1  /  j 

^^  (1   ,     i  =  J 

Then  (by  the  above  assumption)  for  x  7^  0,  J(x)  is  the  Jordan 


tiT' 


form  of  A(x,oo).   Let  1  denote  the  set  of  integers,  s,  for  which 
Xg(x)  =  Xj,  (x) ,  i.e. , 

(1.3)       T  =  s  <=>  X.  (x)  =  X  (x)    for  all  x  in  ^ 

X  3 

It  turns  out  that  no  turning  point  for  x  ^  0   requires 

(l.iia)         fr  (Xj^(x)  -  X  (x))  /  0    for  X  ;^  0 

i.e.  for  X  7^  0  eigenvalues  of  A(x,oo)  which  are  distinct  for 
some  x  ^  0  remain  distinct  for  all  x  7^  0,  and  a  turning  point 
at  X  =  0  requires 

(l.i4.b)         ft'  (X^(x)  -  X,(x))  =  0    for  X  =  0 

i.e.,  at  least  two  eigenvalues  of  A(x,oo)  which  are  distinct  for 
X  7^  0  coalesce  at  x  =  0. 

It  is  well  loiown  that  for  x  ?^  0  a  set  of  n  independent 
formal  solutions  can  be  exhibited  for  (1.1)  restricted  by  (1.2) 
and  (l.l;).   More  precisely,  corresponding  to  the  eigenvalue 
X,  (x)  of  multiplicity  m,  there  are  m,  independent  formal  solu- 
tions each  of  the  form 

00    . 

(1.5)  u(x,p)  =  [exp  ph,  (x)]  511  p"^v,(x) 

^     j=0     J 

where 

(1.6)  hj^(x)  =  I   Xj^(t)dt 

Jo 

One  can  readily  verify  that  the  condition  that  (1.5)  be  a 
formal  solution  to  (l.l)  allows  one  to  determine  the  v.(x)  sue- 
cessively.   (Birkhoff  [1],  Schlesinger  [l],  Noaillon  [1], 


afTi 


Tamarkln  [1]).   Alternately,  one  can  verify  that  the  condition 
that  (1.1)  be  asymptotically  transformed  into  the  diagonal 
equation 

^  w  =  [pXj^(x)  +  . .  .  ]w 

by 

00 
u  = 


_  p"^P:,(x) 
L  j=0     "   . 


w 


allows  one  to  determine  the  P.(x)  successively.   (Hukuhara  [1], 
Turrittin  [2]).   In  any  event,  the  recursion  formula  for  the 
v.(x)  involves  the  pseudo-inverse  of  J(x)  -  X,  (x),  namely 

0  ,   otherwise 


This  pseudo-inverse  and  hence,  as  it  turns  out,  each  v.(x)  (in 
general)  is  singular  at  x  =  0. 

In  order  to  remedy  this  situation  it  is  natural  to  attempt 
to  replace  (1.5)  by  a  more  general  formal  solution  of  the  form 


(1.7)        u(x,p)  =  [exp  ph,  (x)]  2 v.(x,p) 


-22. 

The  fact  that  this  series  is  to  be  an  asymptotic  series  implies 
that  in  an  appropriate  norm,  ||  ||, 


llp^Yvj(x,p)|| 


is  bounded  for  some  x  >  0.  Of  course,  since  x  =  0  is  a  turning 
point  (1.7)  can  represent  the  same  solution  of  (1.1)  in  at  most 
a  sectoral  neighborhood  of  x  =  0.   However,  the  point  of 


introducing  (1.7)  is  to  obtain  a  uniform  expansion  in  the  sense 
that  the  sector  will  include  x  =  0.   (Cherry  [1]). 

Actually,  the  series  in  (1.7)  will  be  constructed  to  be 
convergent  so  that  u  in  (1.7)  will  be  an  actixal  solution  to 
(1.1).   The  evaluation  of  the  series  In  (1.7)  at  x  =  0,  i.e., 

00 

(1.8)  u(0,p)  =  YZ   v.(0,p) 

will  then  be  an  asymptotic  expansion  of  u(0,p)  and  this  consti- 
tutes a  so-called  "connection  formula".   The  step  from  the  con- 
nection to  the  turning  point  to  the  "continuation  around  the 
turning  point"  proceeds  via  Stokes  Multipliers  as  described  in 
Appendix  B. 

In  the  basic  approach  for  obtaining  an  expansion  of  the 
form  in  (1.7)  or,  indeed,  of  the  form  in  (1.5)  one  attempts  a 
perturbation  (or  comparison)  scheme  which  goes  back  to  Liouville 
[1]  and  Steckloff  [l].   This  scheme  regards  the  differential 
equation  (1.1)  as  being  a  perturbation  of  a  comparison  equation. 
The  differential  equation  is  then  converted  into  an  integral 
equation  whose  inhoraogeneous  term  is  a  solution  of  the  comparison 
equation.   ?''ore  precisely,  one  first  defines  a  matrix  G  as  a 
solution  to  a  com.parison  equation 

(1.9)  ^  ^  =  pA(x,p)G 

i.e.,    each  column  of  G  satisfies    the   system 

g^  g  =   pA(x,p)g 

Then  it  is  easily  verified  that  any  solution  of  the  integral 
equation 


(1.10)    u(x,p)  =  G(x,p)c  +  G(x,p)   G"^(t,p)B(t,p)u(t,p)dt 

where  c  is  a  constant  vector  will  also  satisfy  equation  (1.1). 
What  remains  to  be  done  is  to  arrange  matters  so  that  in  an 
appropriate  norm  the  iterant  solution  of  the  integral  equation 

(1.10)  converges  according  to  a  negative  power  of  p.   Por  this 
purpose  define 

(1.11)  ^  E  =  pJE 

(1.12)  V  =  GE"^ 
and  let 

(1.13)  u  =  Vw  exp  ph,  (x) 

Then  it  is  readily  verified  that 

(I.II4.)         /x  ^  =  [p(J  -  \(x))  +  V*"^BV]w 

(Incidentally,  in  the  derivation  of  equation  (1.1[|.)  from 
equations  (1.1),  (1.11),  (1.12),  and  (1.13)  the  diagonal  nature 
of  J  is  irrelevant.) 

The  columns  of  E  must  be  asymptotically  of  the  same  order 
of  magnitude  at  the  connection  point  x  =  0  and  so  one  takes 

(1.15)  E(x,p)  =  exp  p  I   J(t)dt 


V  o 
i.e.. 


:(x,p)  =  <'5^.  exp  ph^(x)i 


as  the  appropriate  solution  of  (1.11). 

In  the  standard  comparison  technique  equation  (l.lLj.)  is 
compared  with 


(1.16)  A.  H  =  p(J  -  Xj^(x))H 

to  obtain  w  as  the  solution  of  the  integral  eauation 

px 


(1.17)   W3  =  6j^3 


where 


+    exp  P(h3j^(x)  -  h3j^(t))  ^  {y-^B^)^M.   dt 

J  x  '^~ 

3 


(1.18)  ^sk^""^  =  ^s^""^  "  \^^^ 
hj^(x)  is  defined  in  (1.6) 

6,  (x)  is  defined  in  (1.2) 
and   the  path  from  x  to  x  is  one  along  which  Re  ph  ,  (x) 

S  3  iC 

decreases. 

The  perturbation  B  is  assumed  to  be  so  small  asymptotically 
that  in  an  appropriate  sector  /(k)  of  the  complex  x-plane 

,-»x 

(1.19)  p^  !   |(V"-'-BV)  .dt|  <  const 

^  X 
s 

for  some  5  >  0.   The  integral  operator  in  (1.17)  is  then  asymp- 
totically contracting  in  (maximum)  norm  since  the  exponential 
terms  in  (1.17)  are  uniformly  bounded  for  all  x  in  Aik)    and  all 
t  on  the  paths.   (The  (maximum)  norm  of  a  scalar  is  a  function 
of  p  and  is  the  maximum  with  respect  to  x  of  its  absolute  value. 
The  norm  of  a  matrix  or  a  vector  is  the  maximum  of  the  norms  of 
its  components.)   The  iterant  solution  of  (1.17)  is  then  major- 
ized in  this  (maximum  norm  by  the  geometric  series 

°°   -5 


const  2. P 

J=0 


n.s 


8 

Suppose  Inequality  (1.19)  is  valid  only  for  6=0  and 
assume,  provisionally,  that 

(1.19)  (V'-'-BV)  .  =  0    when  s  €  k  and  j  e  k 

Then  the  standard  comparison  technique  breaks  down  but  the 
improved  comparison  technique  comes  through  by  the  following 
appraisals. 

A  convergence  factor  is  obtained  for  those  terms  of  (1.17) 
which  contain  exponentials  in  the  integrand  by  an  estimate  of 
the  type  which  is  customary  for  i^aplace  integrals.   This  esti- 
mate is,  for  s  i   k, 


(1.20)      lexp  P(h3j^(x)  -  hgj^(t))(V'^BV)g^.dt|  <  Q.^^'^^ , 


s 


and  it  is  to  be  valid  uniformly  in  xj(k).   A  sufficient  condition 
for  this  estimate  is  that  h,  (x)  have  a  HBlder-continuous  deri- 
vative (which  is  certainly  true  if  A(x,co)  is  analytic)  and 
(V-'-BV)  ,  be  H8lder-integrable,  i.e.,  x(V'-'-3V)  ,  be  HBlder 
continuous.   Then  the  original  procedure  employed  by  Friedrichs 
[1,§22]  is  to  introduce  a  modified  norm  of  the  form 

S  Gk  S  £ik 

(1.21)  iiwiu  =  z:ii^ii  +  P^ziii^ii 

s  s 

where  ||  ||  is  the  above  described  maximum  norm  and  y  is  an 
appropriate  positive  constant.   It  is  then  a  simple  matter  to 
show  that  the  integral  operator  in  (1.17)  is  contracting  in  this 
modified  norm.   A  variant  of  this  final  step,  su.'^gested  by  L. 
^^irenberg,  is  to  rewrite  the  integral  operator  in  (1.17)  for 


-f-,,. 


those  w  with  s  e  k  by  replacing  each  w.  which  occurs  in  the 
s  J 

integrand  by  its  implicit  value  as  given  by  (1.17).   The  integral 
operator  in  the  modified  integral  equation  is  then  contracting 
in  the  original  norm. 

Let  us  now  return  to  (1.19).  Define  projection  operator 
Q  as 

(1.22)  Q  =  \  ( 

0  otherwisej 

and,    for   later  purposes,    the    adjoint  projection 

(1.23)  Q  =  <  ) 

jo  otherwise  j 

Condition   (1.19)    can  now  be   rewritten  as 

(1.21;)  OV'^BVQ  =   0    . 

When   (1.21;)    is   not  valid   instead   of  comparing   (l.lLj.)   with 
(1.16)   we   rather  compare    (l.lij.)   with 

(1.25)  f X  ^  "    [p(J(x)    -   X^ix))    +   QV^BVQJK 

and  it  is  easily  seen  that  the  kernel  of  the  integral  operator 
which  arises  from  this  comparison  will  satisfy  the  decoupling 
condition  implied  by  (1.2i|)   The  inproved  comparison  procedure 
described  above  is  detailed  in  Chapter  II,  and  this  procedure  is 
applied  to  special  cases  in  Chapters  III  and  IV. 

An  implicit  assumption  in  the  preceding  discussion  is  that 
the  connection  formulas  for  the  comparison  equation  are  either 
already  known  or  at  least  are  m.ore  elementary  than  the  connection 


10 

formulas  for  the  original  equation.   In  this  respect  a  natural 
candidate  for  the  role  of  comparison  equation  is 


X   i 


(1.26)  ^G=p^ 

where   each   a,    is    a  non-negative    integer.      By   a  simple    co-ordinate 
transformation   (1.26)    can  be   rewritten  as 

(1.27)  4-  G  =  I      „  •■  G 


dt 


a 
o,  .  t   ,      1  =  n 


A 

The  elements  of  G  can  be  written  e xplicitly  as  (generalized) 

hypergeometric  functions.   (See  appendix  A)   Connection  formula 

for  (1.27)  have  been  determined  by  the  Taylor  series  approach 

only  for  a^  =  0,  1  ^  n  (Scheffe'  [1],  Turrittin  [1])  although  the 

same  approach  would  probably  work  in  the  general  case  of  (1,27). 

In  terms  of  G  from  (1.26)  a  large  class  of  systems  can  be  reduced 

to  a  form  amenable  to  the  improved  comparison  technique.   The 

case  n  =  2  is  relatively  simple  but  rather  important  and  it  is 

discussed  in  detail  in  Chapter  III.   A  typical  result  in  Chapter 

III  is  to  justify  the  conjecture  that  connection  formula  for  the 

scalar  equation 

,2      p 
(1.28)         ^=  [p'^a(x)  +  pb(x)  +  c(x,p)]u 
dx"^ 

where  as  x  ->  0 

a(x)  — >  const  x  ,  a  >  0 
b(x)  ->  const  x^,  P  >  0 


11 


can  be  obtained  in  tern^s  of  connection  formula  for 


(1.29) 

whenever 
(1.30) 


d  g    2   a 
— %•  =  P  X  g 

dx'^ 


a 


P  >  2  -  1 


Equation  (1.26)  can  be  described  as  a  chain  equation  of 
type  zero.   Another  natural  comparison  equation  Is 


f  5 .  J  ,  •  X 
1  1+1, J 


(1.31) 


dx 


C  =  p 


\ 


I  5-  .  X 


n 


5^  n 


i  =  n 


which  can  be  described  as  a  chain  equation  of  type  one.   For  such 
an  equation  the  first  column  of  G  Is  treated  differently  from 
the  other  columns.   This  is  illustrated  in  Chapter  IV  where 
there  is  discussed  a  fourth  order  equat5.on  of  the  type  which 
arises  in  the  theory  of  hydrodynamic  stability. 


II .   Perturbation  of  a  Diagonal  Equation 

It  is  a  surprising  fact  that  a  turning  point  for  a  diffe- 
rential equation  can  be  created  by  a  small  perturbation  of  the 
coefficients.   This  fact  was  demonstrated  by  Friedrichs  [1,§20] 
where,  for  s  implicity ,  the  unperturbed  equation  was 

^1  1     ~  X   0 

=  P 
u^        0   0     u. 


(2.1) 


_d_ 
dx 


U-, 


and,  for  generality,  the  perturbed  equation  was  taken  as 


i) 
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r— 

^' 



P 

X 

0 

+ 

-b^,(x) 

b^2(x)  " 

^l' 

U2^ 

_0 

0  _ 

b23_(x) 

b22(x) 

.^2_ 

(2.2)     I- 
dx 


Friedrichs  [1,§22]  described  a  method  for  determining  inverse 
Stokes  multipliers  for  equation  (2.2)  and  an  elaboration  of  this 
method  forms  the  subject  matter  of  this  chapter.  The  discussion 
here  is  patterned  closely  after  the  discussion  in  Friedrichs 
[1,§21,22]. 

In  analogy  with  equation  (2.1)  the  diagonal  equation  to  be 
perturbed  is 


(2.3) 
where 


^  u(x,p)  =  pJ(x)u(x,p) 


u  = 


'^1 

^2 

« 

• 

• 

u 

n 

_ 

IpI  =  a  large  parameter 


and,  as  in  (1.2),  J(x)  is  a  diagonal  matrix  with  diagonal  entries 
\, ( x) ,  i  =  1,2 , . . . ,n. 
Define 


(2.k) 


'X 

h^(x)  =   X.(t)dt 


The  columns  of  the  matrix 


(2.5) 


G(x,p)  =  (s^  ^  exp  ph.  (x)i 


form  a  set  of  n  independent  solutions  to  (2.3)  and  each  of  these 
solutions  maintains  the  same  exponential  factor  in  a  full 


i      ""Si^l' 
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neighborhood  of  every  point  in  the  complex  x-plane.   Accordingly 
the  unperturbed  equation.  (2.3)  does  not  have  a  turning  point. 
In  general,  the  perturbed  equation 

(2.6)  ^  u(x,p)  =  [pJ(x)  +  P(x,p)]u(x,p) 

will  have  turning  points  even  though  P(x,p)  is  "bounded"  in  its 
dependence  upon  p.   This  fact  will  not  be  demonstrated  here. 
Instead,  a  turning  point  shall  be  presumed  to  occur  at  x  =  0  and 
a  procedure  for  obtaining  a  typical  connection  formula  shall  be 
described. 

Let  the  complex  independent  variable  x  range  over  a  neigh- 
borhood ^   surrounding  x  =  0. 

The  letter  k  shall  be  reserved  to  denote  a  fixed  intep.er 
from  among  the  integers  from  1  to  n  inclusive.   Let  k  denote,  as 
in  (1.3),  the  class  of  integers  s  for  which  \    (x)  =  A-,  (x)  iden- 
tically  and  define  the  projection  operators  Q.  and  Q  as  in  (1.22) 
and  (1.23). 

Let  us  formulate  and  demonstrate  the  key  result  as  a 
theorem  under  convenient  hypotheses  and  subsequently  discuss  the 
reasonableness  of  these  hypotheses.   As  a  preliminary  define 

(2.7)  ^sk^^^  "  ^s^^^  "  ^k^^^ 

Theorem  concerning  the  differential  equation  (2.6). 
Hypotheses 

i)   J(x)  =  \5^  A^(x)| 
ii)   /6(k)  is  a  closed  sector  ( includinp;  x  =  0)  in  /a  with 
the  property  that  there  exists  positive  constants 


Lns 


:!•:  \-. 


Ik 

C  >  0,    C,  ,  and  C   such  that  for  every  s,  /(k)  con- 
tains  points  x  -  x  (k)  and  I'ectifiable  paths  P„(x) 

^  3      S  3 

from  X  to  every  point  x  of  /(k)  such  that  for 
IpI  >  P-  and  for  s  /^  k 

(2.8)  1    lexp  P(hsi,(x)  -  h^^^C  t ) )  |  |Pg^(  t,p)dt  |  <  C^lpl'^^''' 

for  s  e  k 

(2.9)  |Pg^(t,p)dt|  <  Cp 

JX 

s 

where  (T,  C,  ,  and  C   do  not  depend  upon  x. 
iii)   QPQ  =  0,  i.e. 

(2.10)  P   (x,p)  =  0    if  s  e  k  and  r  e  k 

Conclusions 

i)   To  every  k  there  exists  a  solution  of  (2.6),  which 
we  shall  denote  by.a(x,p),  which  is  asymptotic,  uniformly  in 
/(k),  to  column  k  of  the  matrix  G(x,p)  in  (2.5),  i.e., 

max    max    |6,  -  u  (x,p)  exp( -ph,  (x)  )  |  — >  0   as  |p|  -*oo 
s   x  in  l(k)   s^    ^  ^ 

ii)   The  above  solution  u(x,p)  can  be  represented  as  the 
product  of  exp  ph,  (x)  and  an  expansion  which  converges  uniformly 
for  IpI  >  p  ,  p  defined  below  by  (2.26),  and  x  in  Ji{k)    (which 
includes  x  =  0).   This  representation  is 

C0_    /   \ 

(2.11)  u(x,p)  =  [exp  ph,(x)][  +  J__io^    '{x,p)] 

^  r=l 


xioL 
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where  io^{x,p)   =  5^^^ 

and    the  operator  P^  Is  defined  by 


F^u>  =   F(F^"\^)    ,        2   <   r  <  00 
and  for   r  =   1 


s 


(2.12)  (Pv)g  =        [exp   pCh^j^Cx)    -  hgj^(t))]   YL  Pgp(  t ,  p)v^(  t ,  p)dt 

(One   should   note   that    lo  exp   ph,  (x)    is   precisely  column  k   of 
G(x,p))  . 

ili)   In  a  modified  local  norm  the  expansion  (2.11)  is 
asymptotic,  uniformly  for  x  in  i(k),  according  to  integral 
powers  of  p'  ,  i.e., 

(2.13)  max    |  ^   P^.x>  |    <  0,^9'^ '^'■^^'' 
X  in  ^(k)   ^^  "^""^         ^^ 

where 

s  ek  sik 


(2.11+)    Iu(x,p)L  ,  =  YZ   lu  (x,p)|  +  p^'  E:  Iu  (x,p)| 

mod    s    ^  3    ^ 

In  particular,  the  evaluation  of  (2,11)  at  x  =  0, 

00.    _ 
(2.10)         u(0,p)  =  0.  +  ^  (F^(.-0(0,P) 

r=l 

is  an  expansion  which  is  asymptotic  in  the  modified  norm. 
Proof 

Let 
(2.16)  v(x,p)  =  u(x,p)  exp  (-phj^(x))  . 


'I.  ■■■■  ■■: 
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Then  the  Integral  equation  for  v(x,p)  is 


(2.17) 


V  =  i^J  +  Pv 


where   u.^  =  5^^^ 


and    Fv  is  defined  by  (2.12). 

To  present  a  more  detailed  view  of  (2.12)  we  make  use  of 
(2.10)  and  note  that  h  .  (x)  =  0  for  s  e  k.  Accordingly, 

S  xC 


\    TZ   P„^(t,p)v^(t,p)dt  ,    if  s  e  k 


(2.18)   (Fv)   = 


^x  ^ 
s 


sr 


n 


[exp  P(h3i^(x)-hgj^(t))]  XI  Pgj.(t,p)v^(t,p)dt 

s  _ 

'  if  s  /2^  k  . 

Incidentally,  following  the  suggestion  of  L.  Nirenberg, 
the  need  for  a  modified  norm  obviated  if  operator  F  is  modified 


to  be  P  where 


."X 


(Fv)3=)Jx3   - 


I     i_   P^^(t,p)(Fv)^  dt  ,    if  s  e  k 


(Pv)^   of  (2.18)  , 


if  s  ^  k 


The  integral  equation  for  v  is 


A 


V  ^  (V  +   Fv 


but  now  the  iterant  solution  converges  in  the  customary  maximum 
norm. 

Let  us  define  a  global  norm  for  a  scalar  as  follows 
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(2.19)  ||v^(x,p)||  =    max    |v^(x,p)| 

""  X  in  ^(k)   ^ 

Then,  with  the  aid  of  (2.19),  (2.l8)  and  (2.9)  one  readily 
obtains 

(2.20a)        ll(Fv)JI  <  C„  ^  llv^ll    if  s  e  k 

s     p  ^    r 

Similarly,   with  the    aid   of   (2.19),    (2.l8)    and    (2.8)    one   obtains 

(2.20b)  ll(Fv)     II  <  C,p-2*^2i:  llv^ll  if      s^k 

s  ^  r=l        ^ 

Define  a  modified  global  norm  for  a  vector  v  by 

(2.21)  iivii^^^  =  r:iiv3ii  ^-  P^!^iivji 

s  s 

with  Y  as  yet   unspecified.      Clearly   \\  t*>\\      ,    is   finite   for   any 
choice   of  y      Indeed,    ||u.i||      ,    =   1   independent   of  y 

Let  m,    denote    the   multiplicity   of   the   eigenvalue    X,  (x), 
i.e.,   m,    denotes    the   niomber  of  distinct    integers    in  the   class   k. 
Then   (2.21)    is   combined  with   (2.20)    to  give 

r/k  ..   o.<-    n 


11^^ Ld  ^  "^k^p  n:  ii^r"  ^  (--"^k^^hp'^    ^.  "^ 

r  r=i 


which  can  be   rewritten  as 


rek  o^  o  r/k 

llv   , 
r' 


I^^Ld  ^  P""'"^  tc,  zZ  hj  -  (Cip-Y  H-  c^p2^-2Y)pr  ^  II, 


r  r 

where 

(2.22a)  C^   =    (n-mj^)Cj^ 

and 

(2.22b)  C^  =   mj^Cp 


(.!. 
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In  particular,  let  us  select 

(2.23)  Y  =  ^ 

and  define 

{2.2k)  Cp  =  max  (C^p"'^  +  C^,C^) 

then 

(2.25)  llFvlU<Cpp-Mlv|^„^ 

for  p  >  p^. 
o 

Select  p  such  that 
o 

(2.26)  ^pPo'^^  1 

Then  for  p  >  p  the  iterations  of  v  =  ^'  +  Fv  are  dominated 

oo  _^  , 

in   II     II      .,   by  the  convergent   geometric    series  > (C„p~    )  "^ .      In 

mou  •   Q        r    0 

view   of    (2.26)    these    iterations    converge   and  conclusion   (il)    is 
verified.      Since   the   iterations   converge 


N 


-r+l 


II V  -  y-  {F'.j) 


< 


i^^'^'Ao, 


¥^  mod  -         -0- 

F'^o 

But  repeated  application  of  (2.25)  shows  that 

00 

so   that   V   is    asymptotic    to   io -^  >        F  l'-  in  ||     IL    j  •      Moreover, 

max  I II I       ,    <   II  u  IL     , 

.       tf/,  X     '      mod  —  "     "niod 
X  in   >X(k) 

so  that  conclusion  (iii)  is  proved.   Conclusion  (i)  is  a  simple 
consequence  of  (iii),   Q.E.D. 


.•»' 


,   .■■  f      f 


i  y..,.-. 


n:;-i. 


(r'S.S) 
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Better  perspective  on  hypothesis  (iii),  i.e.,  (2.10),  is 
furnished  by  the  following  lenima. 
Lemma  1 .   In  order  that  (2.6)  be  transformed  by 

(2.27)  u  =  Ty 

into  , 

|t  y  =  [pJ(x)  +  P(x,p)]y 


dx 

with 


QPQ  =  0 


it  is  sufficient  that  T  be  defined  by 

(2.28)  T  =  QTQ  +  (5 
and 

(2.29)  ^  ^  =  ^^^-^ 

with  the  non-zero  columns  of  QTQ  independent. 

The  proof  of  this  lemma  proceeds  by  a  direct  verification. 
For  example,  the  existence  of  T~  follows  from  the  definition  of 
T  and  the  fact  that 

(2.30)  T"^JT  =  J 

follows    from   (2.28)    and    the   similar  direct  sum  nature   of  J,    i.e., 

(2.31)  J  =   \Q.  +   QJQ 

Thus  we  see  that  hypothesis  (iii)  is  inessential  and,  in- 
deed, it  can  be  eliminated  by  incorporating  Lemma  1  into  the 
preceding  theorem. 

In  order  to  simplify  the  "connection  to  the  turning"  one 
can  take  the  identity  matrix  as  the  initial  value  for  T  at  x  =  0, 
i.e.. 


■<    ' 


-nl 
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(2.32) 


T  =  Q  +  Q 


at  X  =  0  , 


Of  course  T  still  occurs  in  the  asymptotic  expansion  of  u  and  it 
is  often  desirable  to  seek  the  simplest  T.   For  this  purpose  one 
should  note  that  all  the  elements  of  P  need  not  partake  in  the 
definition  of  T.   Suppose 


(2.33) 


p  =  P°  +  pl 


and   (2.9)    is   replaced   by 


(2.3Ua) 


(2.3l;b) 


|(pl)_dt|    <  C'p-^ 


sr 


|(P°)3rdt|     <C 


and    (2.10)    is    replaced   by 
(2.35) 


QP^Q  =   0    . 


if  s   e  k 


J 


Then  the   preceding   theorem  will   apply   after  obvious  modifications. 
For  example,    the    change   in   (2.2l|.)    is    that  C,    will  be   replaced  by 
C^  +  r^i^C' .      Consequently   it    is   sufficient  to  define  T  by   (2.28) 


and 
(2.36) 


^  T   =   QP^QT 


The  last  item  which  remains  to  be  done  is  to  delimit  a 
sector  i(k)  in  which  (2.8)  is  valid.   This  sector  depends  almost 
entirely  on  the  contours 

Re  ph  .  (x)  =  constant  . 


"•  '.   -^     *        ,v  » 
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The  path  from  x   to  x  is  to  be  one  along  which  Re  ph  ^(x) 
decreases.   In  order  that  x  =  0  be  on  such  a  path,  or,  as  we 
shall  say,  in  order  that  x  =  0  be  accessible  from  x  ,  it  is 


necessary  that  either 


or 


The  contoxirs 


^«  ^k^^s^  ^  ° 


X3=  0 


Re  phgi^(x)  =  0 


are  called  transition  lines  since  as  one  crosses  this  line  the 
relative  dominance  between  exp  ph.  (x)  and  exp  ph  (x)  alters. 
These  transition  lines  issue  from  x  =  0  since  by  (1.6) 

hg(0)  =  h^(0)  =  0  . 

VJhen  h  ,  (x)    has   a   zero   of  order  r    (r  an  integer)    at   x  =   0  then 
the   transition  lines   partition   the   neighborhood   of  x  =   0   into 
2r  so-called  minimum  sectors.      For  a   fixed  x     the   x  which  are 

3 

accessible  forms  a  sector  at  x  =  0  of  at  most  three  adjacent 
minimum  sectors  (Fowler  [l],  Priedrichs  [1,§18]). 

The  estimate  in  (2.8)  can  be  made  as  follows.   Introduce 
the  path  length,  ^,   measured  from  x  as  a  new  variable  of 
integration,  i.e., 

^=  nt)      with   /(x)  =  0  and  /(x  )  =  L 

s 

where  L  is  the  total  path  length.  A  sufficient  condition  for 
(2.8)  is  that  there  exist  constants  k_  with  Re  k,  >  C,  r  >  1, 
kg  >  0,  0  <  1,  such  that  for  all  t  on  the  path  from  x  to  x 
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(2.36)  Re   pi\^U)    -  hgj^(t))    <   -|p|    Re   ]^^i^ 


(2.37)  iPps^^'P)   t^'    ^  M" 


The  estimate   in   (2.8)    is   then  valid  with 

1    -   © 


or  = 


/.QD 


Cj^  =   k2    \        lexp    (-lCj_-('^)||i."^dye| 


(when  L  =  00  it  is  often  useful  to  note  that  in  the  range  where 

A  is  boundedly  distant  from  zero  it  is  sufficient  to  replace  S." 

in  (2.37)  by,  say,  exp  k^£^.) 

In  most  circumstances,  if  a  point  x  is  at  all  accessible 

from  X  then  it  is  accessible  along  a  path  on  which  (2.36)  and 
s 

(2.37)  and,  hence,  (2.8)  is  valid. 

III.   The  Second  Order  Equation 

A.   Introduction 

In  this  chapter  the  turning  point  problem  (see  Chapter  I) 

shall  be  solved  (in  conjunction  with  Chapter  II)  for  certain 

equations  of  the  form 

d  u,     p 
(3.1)  — ~=  [p  a(x)  +  pb(x)  +  c(x,p)]u  . 

dx"^ 

where   p  is  a  large  real  parameter 

X  ranges   over  a    finite  neighborhood   surroiondlng  the 

turning  point 

a(x),  b(x)  and  c(x,p)  are  regular  in  x 

c(x,p)  is  uniformly  bounded. 
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For  convenience,  the  turning  point  shall  be  at  x  =  0.  Denote 

the  order  of  the  zero  of  a(x)  and  b(x)  at  x  =  0  by,  respectively, 

a  and  p. 

In  connection  with  (3.1)  it  is  natural  to  consider  the 

(related)  equation 

2 

(3.2)  ^  S.  =  P^x'^g,  . 

dx 

The  main  point  of  this  chapter  is  to  indicate  a  method  by 
which  the  above  two  equations  can  be  related  whenever 

(3.3)  p  >  I  -  1   . 

The  situation  where  a(x)  is  identically  zero,  i.e.,  a  =  oo, 

will  be  excluded.   This  involves  no  loss  of  generality  for 

l/2 
&=   p  '      can  serve  equally  well  as  the  large  real  parameter. 

In  terms  of  o'  (3*1)  is 

— w-  U-,  =  [o-  b(x)  +  c(x,s-  )]u, 
dx*^  ^  ^ 

and  since  a  turning  point  occurs  at  x  =  0,  b(x)  is  not  identi- 
cally zero. 

The  restriction  (3.3)  is  reasonable  from  the  following 
heuristic  considerations.  The  formal  asymptotic  solutions  of 
(3.1)  and  i3'2)    are,  respectively. 


u,  ^  )[a(x)]  ^  exp  p\  [a(z)]^  dz  I  exp  ^   [a(z)3  ^  b(z)  dj 

^1  1         i    V 


X        1    ^         pX 

J 

g-j^  --v  Ix  J  ^  exp  p  i  [z  J   dz 


rp-^f  'V  .-,->;-.•  -,  _/« 


isn  zl   ^i    (  r,?) 


.^    i:;.'e   O 


i:ox;;p  ;/    ■  ^aJ^itsl  i 


V  -s 


£>fia 


2k 


These  asymptotic  forms  resemble  each  other  down  to  the  turning 
point  provided  [a(x)]"'^'^  b(x)  is  integrable  at  x  =  0  and  (3.3) 
Is  a  sufficient  condition  for  this  integrability . 

An  alternate  heuristic  procedure  is  to  first  rewrite  (3-1) 
in  terms  of  a  stretched  variable  t  (Friedrichs  [2])  where 


to   obtain 


dt 


t   =    p®x      ,        0   >  0 

/    2-20-ae        .a  ^      l-2Q-pO  ,    4-P   x  ^, 

u,    =   (p  at+p  ^     ^o  ...;u^ 


This   resembles    (3.2)   provided  the   term  t      is  dominant  for  all 
possible   stretching,    i.e., 

2   -   20   -   aQ  =   0        implies        1   -   20   -   p0  <  0 

and   this    condition   is   precisely   (3.3).      Thus    (3.3)    may   be   des- 
cribed  as    the   non-transitional   case    in   that   a(x)    is   the   dominant 
term  both  when  x  is  boundedly  distant   from  zero   as  well   as   when 
X   is    infinite simally   close   to   zero. 

For  the  sake   of  clarity  and   for  generality  the   analysis    is 
carried   out    in  terms    of   a  system  of  the    form 


(3.l;a) 
where 


^  u  =    [pA(x)    +  B(x)    +    p"^C(x,p)]u 


u  = 


u. 


u. 


and 


A(x)   = 


0        1 
a(x)      0 


-.'  «t 


■:}tal 


.      -TC 


•M 
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Non-zero  entries  in  C(x,p)  will  be  admissible  provided  that  each 
entry  is  uniformly  bounded.   Similarly,  none  of  the  entries  in 
B(x)  need  be  identically  zero  provided  only  that  each  entry  has 
a  zero  at  x  =  0  of  order  greater   than  or  equal  to  p  vjhere  p  is 
constrained  by  (3  •3). 

The  related  system  for  {3*h^    is  the  matrix  equation 

(3.5a)  ^  ^  =  pA(x)G 

where 

'  0    1 
(3.5b)         AU)  = 

_x    0 

In  order  to  relate  (3.I4-)  and  (3.5)  one  must  first  either 
transform  (3.1+a)  so  that  a(x)  of  (3•i^-b)  becomes  x  or  else  one 
must  transform  (3.5a)  so  that  x  of  (3 •5b)  becomes  a(x).   For 
convenience,  the  former  course  is  selected  here  and  detailed  in 
section  B.   The  system  which  results,  i.e.,  (3.14-)  with  a(x)  re- 
placed by  X  ,  is  called  a  canonical  system. 

In  section  C  we  will  exhibit  solutions  to  a  generalization 
of  (3*5)  which  serve  to  reduce  a  corresponding  generalization  of 
the  above  canonical  system  to  the  form  where  A(x)  is  diagonal, 
B(x)  is  integrable,  and  p~  C(x,p)  is  contracting.   The  further 
discussion  of  this  latter  system  proceeds  as  in  Chapter  II* 

B.   Reduction  to  the  Canonical  System 

This  section  shall  present  in  matrix  form  the  transforma* 
tion  exploited  by  Langer  [1]. 

Let 

X  =  x(x)  ,   u  =  P^v 

o 
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where 


P  = 
o 


0 


0   x' 


'  denotes  differentiation  with  respect  to  x. 

When  u  satisfies  (3.[|.)  the  equation  for  v  is 


(3.6)   x-  ^  V  = 


0      x' 
a(x)    ^ 


T^ 


0   0 


+  P'-'-BP  +  p'-'-p'-^CP 
o   0  "^   o   o 


V 


The  determination  of   x(x)    follows   from  the   requirement   that 


a(x)   =  X   (x    ) 


i.e., 
(3.7) 

where   p  =   a  +  2. 


x(x)   = 


f  j       [a(z)]^  dz 
0  o 


Of  course,  the  implication  is  that  x  =  x(x)  is  a  1-1 
regular  transformation  at  x  =  0.   This  is,  Indeed,  the  case 
since  near  x  =  0,  3!f(x)  approaches 


r^ 

a 

2 

2 

2 
z"^  dz 

"P 

0 

and  this  is  precisely  equal  to  x. 

It  may  happen  that  the  zero  of  x"  is  greater  than  -^  -  1  so 
that  (3«7)  is  already  .n  the  desired  form.   If  this  is  not  the 
case  then  it  turns  out  that  the  x"/x'  term  of  (3.6)  can  be  re- 
moved at  no  expense  to  A(x)  or  B(x)  by 

<|)   0 


(3.8) 


V  = 


-P'-^d*'  4, 


w 


r  .-f ,'      T 


;r.    p.f    ''^    ■^:; 


! 
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where  <j)  x'  =  1.   Consequently,  there  is  no  loss  of  generality  if 
in  (3'k-)    a(x)  is  taken  to  be  x  . 

C .   Diagonalization  of  the  Canonical  System 

In  this  section  we  shall  allow  ourselves  the  luxury  of  dis- 
cussing a  canonical  system  which  is  more  geneiral  than  is  needed 
for  (3.1).   This  system  is 


(3.9) 
where 

(3.10) 


^  u  =  [pA(x)  +  B(x)  +  p'^C(x,p)]u 


A(x)  = 


a. 


0     X 

X  "^   0 


a^  and  a.^   are  positive  integers 
B(x)  has  a  zero  of  order  p  at  x 


=  0 


with 
(3.11) 


p  >  _1  __^_  -  1 


and    C(x,p)  is  uniformly  bounded. 

Incidentally,  under  a  transformation  of  the  form  x  =  2 
where  k  is  a  positive  integer,  a,  ,  a„,  and  (3  change  in  such  a 
manner  that  (3.II)  is  invariant. 

We  first  analyze  the  matrix  equation 


(3.12) 


^G   =  pA(x). 


dx 


For  X  7^  0,  the  Jordan  form  of  A(x)  is  a  diagonal  matrix  which  we 
denote  by  J(x)  and  we  define  E  by 


(3.13) 


:S^  E  =  pjE 
dx     "^ 


;■■-./, 
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Our  goal  is  the  "diagonallzation"  of  (3.9)  according  to 
the  strategy  outlined  in  Chapter  I.   For  this  purpose  it  is  suf- 
ficient to  construct  a  solution  G  to  (3 "IS)  such  that  for  x  in  a 
sector  in  the  x-plane  and  with  V  defined  by 


(3. Ill) 


V  =  GE 


-1 


the  following  three  conditions  are  satisfied. 


/a)   V"-^BV  is  integrable 

b)  p'-^V'-'-CV  is  contracting 

c)  V  is  known  asymptotically  at  x  =  0, 


(3.15) 

With 
(3.16) 

the  equation  for  w  is 


u  =  Vw 


(3.17) 


~  w  =  [pJ  +  V'-'-BV  +  p"-'-V"-^CV]w 


and  when  (3.15)  prevails  the  further  discussion  of  (3.17)  pro- 
ceeds as  in  Chapter  II. 

Condition  (3.l5c)  results  from  the  fact  that  V  can  be 
written  explicitly  in  terms  of  Hankel  functions.   Indeed,  condi- 
tions (3.l5a)  and  (3.l5h)  could  be  demonstrated  from  the  repre- 
sentations of  the  Hankel  functions  as  definite  integrals.   How- 
ever, this  integral  representation  will  be  ignored  so  that  the 
present  discussion  can  be  applied  almost  directly  to  certain 
systems  of  order  n. 

The  parameter  p  is  suppressed  from  (3-12)  and  (3.13)  when 
one  goes  into  the  plane  of  the  stretched  independent  variable  t 
where 
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(3.18) 
Define 


t  =  p^  X 
P 


a,  +  a^  +  2 


G  = 


^2-^1 
.  2p 


<^l-^2 
2p 


.0        P 
The  so-called  stretched  equation  is  then 


(3.20) 


A  G=  A(t)G 


where  A  Is  defined  by  (3.10).   In  terms  of  t  the  solution  of 
(3.13)  is 

0 


(3.21) 


2  I  ri 

E  =  exp  -  t^ 
^  P 


0   -1 


A  line  along  which  Re  t^'   =  0  is  called  a  transition  line 
(see,  for  example,  Friedrichs  [1,§19]).   If  o.^   +  cip  is  not  an 
even  integer  we  have  to  consider  a  doubly  sheeted  t-plane  or 
what  amounts  to  the  same  thing  a  singly  sheeted  t-plane  where 
t  =  t  .   The  variable  x  would  then  correspond  to  any  one  sheet 
of  the  t-plane.   However,  in  order  to  simplify  the  discussion, 
we  shall  suppose  that  the  integer  a  +  a^  is  even.   There  are 
then  p  distinct  transition  lines  and  these  occur  at 


(3.22) 


arg  t 


Tth 


h  =  1 ,3 > . • . >2p-l  . 


2ii 
Consider   the    set   of  p   closed   sectors   each  of   angle    less    than  — 

^  P 

and  each  of  which  is  bisected  by  a  different  transition  line. 
These  sectors  form  an  overlapping  covering  of  the  turning  point. 


•ri>--?- 


;i     \r.X.^) 
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In  this  chapter  the  discussion  shall  be  restricted  to  a  typical 
such  sector  and  this  sector  shall  be  denoted  by  A*      Of  course, 
one  should  keep  in  mind  that  x  is  to  be  in  a  finite  neighborhood 
of  X  =  0  so  that  in  A,    lt|  will  not  grow  more  rapidly  than  p  '^, 
It  is  easily  seen,  see  Appendix  A,  that 

r(l) 


(3.23)   G= 


^ 


+1 


Qp+l 


2t 
2 


r-,,  UK) 
(-^) 


-(^-) 


K 


where  ^  =  —  t^'   and  K  is  a  constant  matrix,  is  a  solution  to 
(3.20) .   (Of  course,  K  will  be  selected  so  as  to  have  a  non- 
vanishing  determinant.)   Define  Cf  by  (3«23)  (with  K  not  yet  spe- 
cified), G  by  (3.19),  E  by  (3.21),  and  V  by  (3.1t|.).   For  each 
sector/^  it  is  a  simple  matter  to  define  K  in  (3.23)  so  as  to 
obtain  the  following  two  estimates  for  V, 


(3.2[^) 


and 


(3.25) 


V  <  const  X 


V   const  p 


^-^2 


°l-°2 


2p 


Let  us  accept  these  estimates  for  the  moment.   In  order  to  com- 
pute V~  we  first  note  that  the  determinant  of  V,  denoted  by 
det  V,  is  a  non-zero  constant  independent  of  x  and  p  since 
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det  V  =  det  G     (by  (3.II4-)  and  (3.21)) 
=  det  G  (by  (3.19)) 

=  constant   (by  (3.20)  since  the  trace  of  A(t)  =  0) 
7^  0   (by  (3.23)  since  the  columns  of  G  are  independent) 

Consequently,  the  elements  of  V"   are  a  permutation  of  the 
elements  of  V  so  that  estimates  analogous  to  (3.21+)  and  (3.25) 
are  valid  for  V~   and  finally 


(3.26) 

and 

(3.27) 


V'-^BV  <  const  X 


P- 


^1-^2 


1^1-^2 


p'-'-V'-'-CV  <  const  p'   P 


This  already  verifies  (3-15)  since  by  (3.11) 


o^-a^ 


and  by  (3.l8) 


^-^2 


>  -1 


-  1  <  0 


The  only  item  which  remains  to  be  completed  is  the  deriva- 
tion of  the  estimates  (3.2i|)  and  (3.25). 

The  behavior  of  G  in  a  finite  neighborhood  of  t  =  0  is 
already  described  by  (3.20).   In  order  to  estimate  the  behavior 
of  G  outside  a  finite  neighborhood  of  t  =  0  let 


(3.28) 


G  = 


a^-a^ 


0 


^2-°l 
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and 
(3.29) 


Y  = 


1   1 
1   -1 


T 


Corresponding  to  (3.20)  for  G,  the  equation  for  Y  is 


dt 


Y  = 


Vz 

2 

'o 

1 

^1-^2 

i 

0 

1 

0_ 

i^t 

jD 

-1 

and  the  equation  for  T  is 


(3.30) 


d       m 

^      2 

1 

o" 

[        °l-''2 

0 

1 

dt     ^     = 

t 

0 

-1_ 

■       5t- 

1 

0 

From  well  known  theory  (see,  for  example,  Priedrichs 
[1,§13])  (3.30)  has  a  solution  T  whose  asymptotic  form  for  t  in 
^  but  t  large  is  a  constant  multiplied  by  the  matrix  E  of  (3*21). 
Select  K  in  (3.23)  so  that  this  T  is  the  solution  of  (3.30)  which 
corresponds  to  (3.23).   One  should  note  that  K  can  be  determined 
explicitly  since  the  continuation  problem  for  the  Bessel  equation 
has  already  been  solved.   (Watson  [l,Chap.7]). 

From  the  above  selection  of  K  and  by  (3.28)  and  (3.29) 


(3.31a)    GE'-*-  <  const 


.~ir" 


a^.a^ 


0 


t~r 


1    1 
1   -1 


,   for  It  I  >  t^ 


and  by  (3.20) 
(3.31b) 


^E"  <  const  ,   for  |t|  <  t 


where  t   is  a  fixed,  but  sufficiently  large,  real  number. 


'  o' 
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For  definltene33,    let    °-p-\   >  0«      Then,    from   (3.31) 


(3.32)  GE'-*-   <   const 


or,    since   by    (3.l8) 


1        0 


"2-^1 


0        1+t^f 


1        1 
1      -1 


,      all  t   in  /t 


(3.33) 


p  ^       t       ^  =    X       ^ 


(3.19)    I3  combined  with   (3.32)   to  form 


(3.31;) 


GE"      <  const 


p  '^ 

0 

0 

1. 

1      1 
1    -1 


all  t    in   i. 


which,  with  (3.11;),  results  in  {3-2$). 
On  the  other  hand,  from  (3.31) 


(3.35) 


GE   <  const 


°l-^2 


0 


1+t 


[i   i" 

Cp-a^  [1      -1_ 
-Tr- 


ill t  in  i 


and  this  is  combined  with  (3.33)  and  (3.19)  to  for 


m 


(3.36) 


GE"   <  const 


a^-a^ 

X      ^ 

0 

i 

1' 

P 

1_ 

1 

-i_ 

which,  with  (3.114-),  results  in  (3.2i|.). 


<■  r 


J' 


3k 


IV.   A  Fourth  Order  Equation 


In  this  chapter  we  shall  indicate  how  one  can  obtain  con- 
nection formula  for  an  equation  of  the  form 

{k.D      u^^^^  -  p^([a(x)  +  p-H^(x)+p'^c^(x,p)]u^^^^ 

I  o         o 

+  [b^(x)+p"^c^(x,p)]u^^^  +  (b2(x)  +  p'-^C2(x,p)]uj  =  0 
in  terms  of  connection  formula  for 


(U.2) 


g^i-)  -  p2a(x)g(-  )  =  0 


when  at  x  =  0,  a(x)  has  precisely  a  simple  zero  while  each  of 
the  terms  b  (x),  b,(x),  and  bp(x)  has  at  least  a  simple  zero. 
The  terms  c  (x,p),  c, (x,p),  and  C2(x,p)  represent  terms  which 
are  bounded  uniformly  in  x  and  p.   Equation  (i]-.l)  differs  from 
a  classical  equation  which  arises  in  the  theory  of  hydrodynamic 
stability  in  the  above  requirement  that  bp(0)  =  0  and  (I4..2)  is 
correspondingly  simpler  than  the  customary  related  equation 
(Wasow  [1],  Lin  [l]). 

Equation  ik-'X)    is  the  first  component  equation  of  the 
system 


(1*.2)   ^u  = 


0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

+ 

0 

0 

0 

a{x) 

0 

b 

1 

0 
0 


0 
0 
0 


|b2(x)   b^(x)   bQ(x)   0 


+  P 


-1 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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p-^c. 

^0 

0 

u 


■'■■'i  "xs-tnO  II 


:■.:•  I 


lo'i 


ejux\ 


i-  ^ 


«*.   V   •    J-  " 


■-ioiij.:;.. 


(Q.?r)    '^ 


i.j- 11.  ^-i.^ 


yr,x> 


3S 


It    is    a   routine  matter   to   reduce    (l+.l)    or   (1;.2)    to   an 
equation  of   the   same   form  but  where    a(x)    is   replaced  by   x.      (See 
for  example,   Chapter   III,§B.)      For  simplicity,   this   will  be 
assumed  to  have  been  done  and  we  will   start  with  the   equation 


(U.3)       ^u  = 


0  0      0      0 

0  0      10 

0  0      0      1 

0  0x0 


+  B(x) 


+   p"^C(x,p) 


u 


None  of  the  entries  in  B(x)  or  C(x,p)  need  be  identically 
zero  for  the  subsequent  analysis.   The  only  restriction  on  B(x) 
is  that  the  entries  in  the  last  row  vanish  at  x  =  0.  Each  entry 
in  C(x,p)  is  required  to  be  bounded  uniformly  in  x  and  p  and, 
moreover,  p  '^   times  each  entry  in  the  second  column  of  C(x,p) 
is  to  be  uniformly  bounded. 

The  related  system  for  ([j-.3)  is  the  matrix  equation 


^h.k) 


0 

0 

0 

0 

dx  " 

P 

0 
0 

0 
0 

1 

0 

0 

1 

0 

0 

X 

0 

H 


and  this,  in  turn  is  related  in  an  obvious  manner  to  the  system 


namely, 
(I+.6) 


dx 


H  = 


0 

1 

0 

p 

0 

0 

1 

Si 

X 

0_ 

1 

0 

0 

_0  ' 

0 

0 

G 

0 

G 
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Prom  (I4..5)  one  obtains  the  parameter  free  stretched  system 


(I;. 7)                      ^G  = 

0 
0 

1 

0 

0 

1 

G 

_ 

.0 

t 

0  _ 

by  the  definitions 

1 

0 

1 

0 

ik'Q)                                    G  = 

0 

p^ 

0 

_0 

0 

1 

G 


and 
(I;. 9) 


t  =  p2/3... 


The  solution  to  (i+.T)  can  be  written  explicitly  in  terms 
of  Airy  functions  Ai(t)  and  Bi(t),  namely, 


(U.IO) 


G  = 


p  '^   i  Ai(z)dz 


•'O 


.t 


JO 


0 


Ai(t) 


Bi(z)dz 


■Bi(t) 


0    I^Ai(t) 


--^Bi(t) 


0 
0 


0 


K 


where  K  is  a  constant  matrix  (not  yet  specified). 

It  is  clear  from  (1|.7)  that  G  of  ([j..lO)  is  bounded  in  a 
finite  neighborhood  of  t  =  0.   In  order  to  estimate  G  for  t 
large  let 


(U.ll) 


t 

1 

0 

0 

^  = 

0 

1 

t-5 

0 

0 

0 

1 

t5 

.1  -  ^  :t 


.1') 


L 
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The  equation  for  Y  is  then 


dt 


Y  = 


1 
.2 


0  10 
0  0  1 
0    10 


ht 


0 

o" 

1 

0 

0 

-1. 

Finally,  let 

(1+.12) 

to  obtain 


(1,.13)    ^T 


Y  = 


1 
2 


1 
0 
0 

0 

1 

0 


1  1 
1  -1 
1   1 


8t 


5 

0 
-2 


5' 

-2 

0_ 

The  first  column  of  T  is  selected  as 


1  3 


1 

0 
0 


so    as    to  be  consistent   with   (I4..IO).      Define   E   as 


(k-lk) 


E  =  |t2 


0  0  0 
0  10 
0        0-1 


A  transition  line  is  defined  as  a  line  along  which 
Re  t  '   =  0.   Let  ^  be  a  closed  sector  which  contains  one  and 
only  one  transition  line. 

For  t  large  (i]..13)  has  a  solution  T  whose  last  two  columns 
are  asymptotic  to  a  constant  multiplied  by  the  last  two  columns 
of  E.   Select  i^   in  (i|.10)  so  that  this  T  is  the  solution  to 
(1;.13)  which  corresponds  to  (I4..IO). 

Let  [QE'  ]  denote  the  last  two  columns  of  GE'  .   Then  in  A 


:&:■ 
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{i4..l5a)   [GE"-'-]  <  const 


i 


0     0 

1 


0 


t"*- 


t  < 

o  ~ 


ik'lSh) 


[GE"-^]  <  const   ,    ltl  <  t 


where  t   is  a  fixed,  but  sufficiently  large,  real  number.   The 
first  global  estimate  for  GE~  which  follows  from  (i|.l5)  is 


(i4..l6)   [GE'-*-]  <  const 


1 


0 


0 


0 

1 

^   0 


0 


1 

t^ 


1  1 
1  -1 
1   1 


all  t  in 


i.e.,  from  (i;.8)  and  ([|..9) 


(l|.17)   [GE"-'-]  <  const 


3 

0 

0 

0 

1 

0 

0 

0 

1 

x5. 

1 

1  1 
1  -1 
1  1 


,  all  t 

in  J 


The  second  global  estimate  which  follows  from  (Li-.lS)  is 


(i+.l8) 


[GE'-'-]  <  const 


0 
1 


0 
0 


0    0    t 


1 

-1 
1 


all  t  in  y^ 


i.e.,  from  (l+.S)  and  (li..9) 
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(14..19)      [GE"-'-]    <  const 


1 
2 


0 


0 


1 

Z 


0 

] 

p' 

0 


1       1 

IT      'Z 
t  +  p   J 


,      all   t    in  // 


In  order   to  carry   over  these    results   to  H  of   (i^-.ij.)   define 


E  = 


where  E  is  defined  by  (I4..IL1.),  and  define 


1 

0 

0 

0 

0  ■ 

0 

E 

0 

(1;.21) 


V  =  HE"-*-  . 


Then  iU-.ll)    results  in  the  estimates 


(i+.22) 


V  <  const 


10   0 


0   1 

0    0 


while  (I4..I9)  results  in 


(U.23) 


V  <  const 


0 


0   0 


0    0 


X 


X 


3 


1 


0   0   l+x' 


1 


0    0 


1 

1 

z 


•X 


3 

1 


1+x'- 


1 


1 
Z 


-P 


1   1 


1      1 

z    u 

+  x^ 


,-1 


In  order  to  estimate  V  one  first  notes  that  the  determi- 
nant of  V,  denoted  by  det  V,  is  a  non  zero  constant  (independent 
of  X  and  p)  since 


'  ■  (-*>• '  •  '"I    '■ 


ko 


det  V  =  det  H 


(by  {It-.li].)  and  ([|..21)) 

(by  (i;.6)) 

p-'  -  det  (i        (by  (l|.8) ) 

=  (det  K)   Wronskian  of  (Aj^(  t )  ,B^(  t)  )  (by  (U-IC)) 
=  constant  j/   0  , 


=  det  G 
-  .1/2 


det  G 


Consequently,  the  estimate  for  V~   is  the  same  as  the 
estimate  for  the  matrix  which  is  obtained  from  V  by  first  repla- 
cing each  element  by  its  cof actor  and  then  interchanging  rows 
and  columns.   Thus,  from  (I4..22),  one  estimates 


ik'^h) 


V~  <  const 


while  from  ik'23) 


ik'2S) 


Define 


(1+.26) 


V"  <  const 


J  = 


0 


0 
0 


0 


0 
1 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 


0 

1 

1+x 

] 


-X 


1 


.-1 

1 

1 


0 


-p 


1   1 

+  x^ 


1 
Z 


0 
0 


1 


0 


0 
0 

0 
1 


-p 


1 

6 


-  X 


1 

z 


Kiifi-^,.   . 


^1'- 


and 

(l|.27)  w  =  Vu 

where  V   is   defined  by   (I|..21).      Then 

(i|..28)  ^  ^  ^    ^P*^  "^  V'-'-BV   +   p"-^V"-'-CV]u 

In  the  following  we  restrict  ourselves  to  a  finite  neigh- 
borhood of  X  =  0.   Prom  (I4..2I4.)  and  the  fact  that  the  last  row  of 
B  has  a  zero  at  x  =  0  It  follows  that  V"  B  is  bounded.   This 
fact,  in  combination  with  il^.,22)    shows  that 

_3 
(i]..29)  V'^BV  <  X  ^ 

which  is  certainly  integrable.   Finally  from  (I4.23)  and  the  fact 
that  p  '-^  times  the  second  column  of  C  is  bounded  it  follows  that 

1 

CV  <  p^ 

This  fact  in  combination  with  (I4..25)  shoxjs  that 

(I4-.30)  p'-'-v'-^cv  <  p  ^ 

which  is  certainly  contracting. 

The  estimates  (1+.29)  and  (14..30)  together  with  the  obvious 
fact  that  J  of  (I4..26)  is  diagonal  show  that  the  further  discus- 
sion of  (l|.28)  proceeds  as  in  Chapter  II. 
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V.   Conclusions 

The  classical  work  on  the  turning  point  problem  has  relied 
upon  a  canonical  system  which  is  a  perturbed  diagonal  system 
where  the  integral  of  the  perturbation  term  is  asymptotically 
small.   A  more  appropriate  foundation  has  been  elaborated  in 
this  report,  namely  a  canonical  system  which  is  a  perturbed  dia- 
gonal system  where  the  inte?^ral  of  the  perturbation  term  is  mere- 
ly bounded.   In  order  to  clarify  the  scope  of  this  latter  cano- 
nical system  it  is  important  to  delimit  the  class  of  perturba- 
tions which  are  "accessible"  with  respect  to  a  given  comparison 
system.   For  example,  if  the  comparison  system  is 

^G(x,p)  =  pA(x,p)G(x,p) 

then  the  fact  that  the  system 

^  u(x,p)  =  [pA(x,p)  +  B(x,p)]u(x,p) 

is  reducible  to  the  improved  canonical  system  via  solutions  of 
the  comparison  system  will  imply  a  condition  on  B(x,oo  )  and  a 
condition  on  B(x,p)-B(x,oo  ) .  The  precise  nature  of  these  condi- 
tions have  been  indicated  (in  Chapters  III  and  IV)  when  the  com- 
parison equation  has  the  form  of  certain  chain  systems.  The 
results  in  Chapter  III  on  a  second  order  equation  merits  special 
attention  since  it  provides  a  good  perspective  on  the  many 
papers  previously  published  on  this  problem  (See,  for  example, 
Langer  [3],  Jeffreys  [1]). 
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Appendix  A;   The  (General)  Hypergeometric  Equation 

An  equation  of  the  form 

n  ^11^ 


(A-1)        TT  (©  +  P.  )y  =  Yt*"  Vi   (9  +  a  )y 
1=1       ^         1=1       ^ 


where 


Q   =  ^ 
^   dt 


Is  denoted  as  a  (general)  hype rgeome trie  equation. 

Related  to  this  equation  are  the  (general)  hypergeometric 
series  which  are  defined  in  the  following  notation,  [Watson, 
p.  100]. 

00  (a;L)n---(%)n  t^ 


(A-2)    pPq(a^,...,ap;p^,...,Pq;t)  =  ZI  X^)^.  .  .  (  0^)^  ET 

where  (a)   =  a( a+1) . . . ( an-1) ;  (a)^  =  1.   For  this  series  the 
following  lemma  is  readily  verified. 

Lemma  1:   F  (a  , ...,a  ;P  , ...,p  ;t)  satisfies  the  equation 

(A-3)         G(0+p-l)...(Q+p  -l)y  =  t(Q+a^)...(Q+ap)y 

Lemma  1  does  not  yet  apply  to  (A-1)  and  In  order  to  remedy  this 
we  state  two  more  trivial  lemmas: 

Lemma  2 ;   f(0)t''y  =  t''f(0+j)y  where  f(0)  is  a  polynomial  with 
constant  coefficients. 

Lemma  3:   If  z  =  t*^  then  k.Q^  =  Q^   where  0^  =  t  ^  and  ©^  ~  ^  37' 
By  means  of  these  last  two  lemmas  we  now  convert  (A-1)  into 
(A-3)  as  follows.   First  y  of  (A-1)  is  transformed  according  to 

(A-i;)  y  =  t  -J  V 


"^.O    R' 


.;  ;•■ 


kk 

From  Lemma  2  the  resulting  equation  for  v  is 


n  .,  w 

(A-5)      TT  (©  -  p.  +  P.)v  =  Yt'  TT  (0  +  a  +  p  )v 


Next,  with  the  aid  of  Lemma  3,  (A-5)  is  rewritten  in  terms  of 
(A-6)  z  =  ^^^— 


^n-m 


and  the  new  equation  is 


IT  (9„  +  -~~^)y  =  z  TT  (©  +  -~-^) 


(A-7)      TT  (9,  +  -~~^)y 


1=1   z     K    -     i^i   z 


Finally,  we  define  p.  for  i  >  n 


as 


Pi+n  =  Pi 
so  that  equation  (A-7)  can  be  rewritten  as 

n-1       -Pi+1+P^       40.        V^i 

(A-8)   Q^  TT  (e^  +  — ^-^)y  =  2  TT  (92  +  -^)y 

This    equation   is   precisely   in   the    form  covered  by  Lem.^a  1    so 
that    one  can  exhibit    the    following  n  solutions   to   (A-1). 

(A-9)  yj  =  t   Vn.l<°lJ'---^TirPlj"--Pmj'^^ 


a,+p.  -Pi+i+Pi 

i^  -  ^^     arid      Pij  =        •'    i;     -^  +   1,      j=l,2,...,n. 


where   a,  .  =  — r—^     and   p 


These  solutions  will  be  distinct  provided  that  none  of  the 
quantities  p./K.  differ  by  zero  or  an  integer.   Moreover,  (A-9) 
will  represent  actual  solutions  (and  not  just  formal  solutions) 
provided  n  >  m.  When  n  >  m  then  the  series  converge  for  all  t. 


\'    WSJ 


^  -». 


x-r. 


One   should   note   that    the   above   results   give   a  simple   suffi- 
cient  condition  in  order  that  n   independent   regular  solutions 
exist   at    the   s ingular  point   t   =   0.      This   sufficient    condition  is 
that  n   >  m,    p.    be   non-nepative   integers,    and   none   of   tlie    p. 
differ  by   a  multiple   of    K. 

Let  us   now   show   the  hype rgeome trie   character   of  the 
equation, 

(A-10)  fr  <4-  fe)y  =  y 


i  .e. 


t 

n     -d+ti.  ) 
TT  (t     ^  e)y  =  y 
i=l 

By  repeated  application  of  the  shifting  formula  of  Lemma  2  one 
can  rewrite  this  equation  as 

(A-11)  TT  (e  -  «S'^Jy  =  t''''  y 


'►C  =  211  ( 1  ■*"  l^i )  ' 


where  (S-^  =  >   (1  +  ix.  )  ,    o„  =  1 

i=l       ^ 

Clearly  this  is  a  (general)  hypergeometric  equation.   Moreover, 
from  the  criterion  discussed  above,  (A-10)  will  have  n  indepen- 
dent regular  solutions  at  t  =  0  provided  the  \i.    are  distinct 
non-negative  integers. 

When  n  =  2  (A-10)  or,  more  generally  the  system 

d       ,\ 

dt  Si  -  t   g2 

(A-12)  ^ 

^       a        -  t-  2  „ 

dt  ^2  -  t   Si 


il'i 


-ilV. 


1X9 


t  rsB'-i 


l|6 


can  be   solved   explicitly   in  terms   of  Bessel   fionctions    J^^(t) 
where    [Watson,    p. 83] 


(A-13) 

and,  consequently, 
(A-lli.) 


(e  +  v)j^(t)  =  tj^_^(t) 


(e  -  v)j^(t)  =  -tJ^^;^^*) 


(©2  -  v2)j^^(t)  =  -t^j^^(t) 


with  the  aid  of  Lemma  2,  (A-I3)  can  be  rewritten  as 


(A-15) 


et^j^^(t)  =  ±t^-'j^(,.i)(t) 


(A-16) 
where 

Let 

(A-17) 

and 
(A-18) 


The  equation  for  g, ,  corresponding  to  (A-11),  is 
Q(9  -  (a^  +  l))g^  =  tP  g^ 

p  =  a,  +  a2  +  2 


a^+1 


gi  =  t 


h. 


z  =  t 


2 


then,  from  (A-I6)  and  Lemmas  2  and  3 


e 


a  +1  2 
P 


4  =  -<^-)'  '^1 


so  that  by  (A-ll(.) 
I    (A-19) 


hi 


In  other  words, 


a^+1 


(A-20) 


Si 


=  z 


a  +1  ^  p 
+  -i— ) 
-   P 


\~:r   z) 


For  the  determination  of  go*  write  the  first  equation  of 
(A-12)  as 


(A-21) 


g2   2  ^ 


^zSl 


and  apply  (A-l5)  to  obtain 


a^.l 


(A-22) 


S2 


+  £  z 
-  2  ^ 


^2i   V 

+(-4--) 


Appendix  B:   Stokes  Multipliers 


Let  U  denote  a  matrix  each  of  whose  columns  satisfies  (1.1), 


i.e . , 
(B-1) 


dx 


U  =  [pA(x,p)  +  B(x,p)]U 


Let  E(x,p)  be  defined  as  in  (l.l^),  i.e.. 


(B-2) 


E(x,p)  =  exp  p 


J(t)dt 


0 


where    J(x)    is    a  diagonal   matrix  vjhich  for  x   in   /?  but   x  7<^  0   is 
the    Jordan  form  of  A(x,oo).      Then   (B-1)    has    a   formal  solution, 
for   X  7^  0,    of   the    form 


•Hi 


(B-3)  U  = 


I4O 
00 


p"^V.(x) 


E(x,p) 


where  the  columns  of  V  (x)  are  independent.   Let  (^denote  an 
annulus  of  f.   boundedly  d5.stant  from  the  turning  point  x  =  0. 
Then  ^-  can  be  covered  by  overlapping  sectors  i.    in  each  of  which 
there  exists  an  actual  solution  of  (B-l)  which  is  asymptotic  to 

the  formal  solution  uniformly  in  &..      Denote  this  actual  solution 

i  ^ 

by  U   or,  where  convenient,  by  U   with  the  superscript  i  under- 
stood.  Let  U  denote  an  actual  solution  of  (B-l)  which  reduces 
o 

to  the  identity  matrix  at  x  =  0.   Prom  the  linearity  of  (B-l) 
and  the  linear  independence  of  the  columns  of  U_^  and  U  it 

^  OD        o 

follows   that  any  solution  U  of   (B-l)   can  be   represented  as 

(B-k)  U  =   U  M    (U)    =   U     M      (U) 

^  00'  00    00       ' 

where  the  matrices  M  (U)  and  M   (U)  are  independent  of  x.   The 
elements  of  M   (U)  are  called  Stokes  multipliers  of  U  with  res- 
pect to  U   while  the  elements  of  M  (U)  are  called  inverse 
'^        00  o   ■ 

Stokes  multipliers  of  U  with  respect  to  U   (Turrittin  [1]). 

One  first  evaluates  inverse  Stokes  multipliers  of  U_^  as 

00 

follows.      Apply    (B-l|)    to   U        and   evaluate    each  term  at   x  =   0   to 
obtain 

(B-5)  M^(U^)    =   U^      at   x=   0    , 

Thus  the  columns  of  M„(U   )  are  connection  formula  and  a  typical 
column  is  evaluated  iri  the  manner  indicated  by  the  theorem  in 
Chapter  II. 

Incidentally,  Stokes  multipliers  for  U  can  be  determined 
from  inverse  Stokes  for  U   as  follows.   Apply  (B-l;)  to  U  and 


'_  J  "Ji>."t"    O":  ;"■    ' 


1^9 

evaluate   each  term  at   x  =   0   to  obtain,   with   the    aid   of   (B-5), 

(B-6)  M      (U   )    =    [M   (U      )]"^ 

'  coo  o^oo 

Finally,  one  notes  from  (B-i|)  and  (B-5) 

(B-7)  M  (IT   )N   (U)  =  M  (U   )M   (U) 

'  ooo'co'    ooo'co' 

Thus,  from  a  knowledge  of  Stokes  multipliers  for  i.  one  obtains 

/•J 
Stokes  multipliers    for  ^ ,   as 

j  i 

(B-8)  ^00^"^    "  ^  ^00^^^ 

where 

K  =    [M    (U      )]    -^  M   (U      ) 

O        00  O        00 


"s    ?•,..  oi:S,>  ■/■ 
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